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1.0 


IHTRODUCTION 


Key  parameters  in  the  assessment  of  effectiveness  of  high  energy  laser 
systems  arc  the  peak  intensity  and  spot  size  in  the  focal  plane.  These  tvro 
features  arc  strongly  influenced  by  interactions  between  the  laser  beam  and 
the  atmosphere  in  the  propagation  path.  Die  primary  interaction  of  concern 
here  is  heating  of  the  atmosphere  due  to  absorption  and  the  hydrodynamics  of 
the  heating  induced  density  changes  in  the  air.  I^ecification  of  parameters 
such  as  peak  intensity  then  becomes  a problem  of  specifying  in  detail  the 
nature  of  this  self -induced  atmospheric  lens  through  vrtiich  the  beam  must  pass. 

In  addition  to  molecular  absorption,  which  sets  the  scale  on  how  much 
energy  is  available  to  cause  distortions  in  the  medium,  the  other  critical 
factor  is  the  rate  at  which  air  flows  across  the  beam.  Die  amount  of  heating 
experienced  by  a parcel  of  air,  sutid  hence  its  distorting  power  as  a lens,  is 
determined  by  the  amount  of  time  it  resides  within  the  laser  beam.  The  extent 
to  which  the  laser  beam  itself  is  distorted  is  then  a strong  function  of  the 
details  of  the  flow  field.  Important  limiting  assunptions  have  been  made  in 
theoreticail  analyses  of  this  effect  which  is  commonly  called  thermal  blooming. 
These  generally  include: 

• Steady  state  conditions  prevail 

• Flow  rates  are  constant  (space*  and  time) 

• Time  dependence  only  for  zero  flow 

Actual  high  energy  laser  engagements  will  generally  violate  one  or  more  of 
these  assumptions.  The  hydrodynami cal  formulations  which  describe  the  laser  - 

« 

lb  account  for  slew,  the  flow  velocity  is  normally  allowed  to  increase 
linearly  with  distance  along  the  beam. 


atmo;;])ii(.ri'  interactions  must,  for  theao  cncaccments,  be  cast  in  a form  such 


tliat  the  a33uji5)tions  listed  above  are  not  required.  The  study  reported  here 
dealt  with  these  limitations  and  consists  of  two  distinct  segments,  namely 

• Analytical  studies  of  hydrodynamics 

• Wave  optics  propagation  calculations 

The  hydrodynamical  formulation  chosen  must  be  general  enough  to  permit  des- 
cription of  cases  which  include 

• Non  coplanar  slew  and  wind 

• Stagnation  zones 

• Accelerated  beams 

• Variable  winds  along  the  beam 

all  of  which  are  likely  to  be  encountered  in  typical  engagements.  Given 
hydrodynamical  solutions  applicable  to  air  flow  conditions  as  above,  the 
effects  on  laser  beam  properties  (e.g.  focal  plane  irradiance)  must  be 
determined  by  performing  wave  optics  propagation  calculations  for  a beam 
transmitted  throiigh  density  profiles  as  determined  self  consistently  for  the 
local  laser  beam-air  interaction. 

In  Section  2.0,  we  show  the  detailed  solution  to  the  hydrodynamical 
equations  for  one  spatial  dimension.  In  this  formulation,  the  wind  velocity 
may  vary  in  magnitude  with  time  and  density  perturbations  in  the  medium  may 
be  tracked  from  beam  turn  on  to  the  stea^  state  - if,  in  fact,  the  steady 
state  is  over  attained  in  times  of  interest.  The  primary  advantage  of  this 
derivation  is  its  physically  intuitive  foundation.  The  various  factors  which 
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drive  the  blooming  process  separate  into  couponents  which  have  clear  physical 
interpretations.  In  the  multi -dimensional  theory,  given  in  Section  5.0,  this 
ability  to  clearly  identify  specific  physical  interactions  is  frequently  lost 
in  the  details  of  the  mathematics.  Bie  formulation,  however,  permits  flow 
velocities  across  the  beam  to  be  con^sletely  arbitrary  in  magnitude,  direction, 
and  time  dependence.  Con^jutations  of  density  contours  and  focal  plane 
irradianco  profiles,  for  this  more  detailed  multi -dimansioiml  theory,  ajre 
shown  in  Section  4.0. 
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2.0 


SOmTION  fyDR  ONE  SPATIAL  DIMENSION 


For  predicting  thermal  blocsning  of  a laser  beam  propagating  in  the 
atmosphere,  it  is  usually  assumed  that  the  unperturbed  airflow  through  the 
beam  is  Independent  of  time.  This  assumption  is  not  valid  for  a slewed  laser 
beam  in  the  presence  of  wind,  as  well  as  in  the  obvious  cases  of  variable 
wind  velocity  or  slew  rate. 

Considerable  theoretical  work  has  been  done  on  the  problem  of  time- 
dependent  changes  in  the  density  of  a medium  through  which  a constant-velocity 
laser  beam  is  propagating.  Recently,  EUinwood  and  Mirels^^^  calculated 
two-dimensional  density  profiles  within  a slewed  laser  beam  for  both  subsonic 
and  supersonic  wind  transverse  to  the  beam,  using  linearized  hydrodynamic 
eqxiations.  Hayes has  addressed  the  non-linear  effects  to  be  expected 

when  the  slew  rate  is  near  the  sound  speed,  using  one-dimensional  equations. 

(4) 

Fleck  et  al.  ' have  developed  an  elaborate  "four-dimensional"  computer 
code  for  calculating  thermal  blooming,  based  essentially  on  steady-state  equa- 
tions for  the  density,  but  capable  of  treating  supersonic  as  well  as  subsonic 
wind  speeds.  All  these  treatments  assume  both  the  wind  speed  and  direction 
to  be  constant  at  a given  position  along  the  beam. 

We  here  present  the  general  time-dependent  solution  to  the  linearized 
one-dimensional  hydrodynamic  equations  for  a beam  with  an  arbitrary,  time- 
dependent  intensity  profile  auid  which  moves  with  variable  velocity.  For 
large-diameter  beams  and  relatively  short  irradiation  times,  the  effects  of 


conductivity  and  convection  are  expected  to  be  small  and  have  been  neglected 
(cf.  Rr  f.  ‘j) . Particular  solutions  for  accelerated,  Gaussian  CW  beauns 
au'e  compared  with  our  analytic  solution  for  a Gaussian  beam  moving  at  constant 
velocity.  Greatest  differences  between  the  solutions  occur  for  beams  which  are 
moving  at  mach  numbers  near  zero  or  unity.  Density  changes  within  an  accelerated 
beam  can  be  either  larger  or  smaller  than  for  an  unaccelerated  beam  moving  at 
the  same  speed. 


For  the  purpose  of  obtaining  the  density  perturbation  we  assume  the 
unperturbed  medium  to  be  at  rest  while  the  laser  beam  moves  through  it. 

Of  course,  this  Is  fully  equivalent  to  the  case  of  wind  blovang  past  a stationary 
beam,  and  is  a good  approximation  for  a slewed  beam  as  long  as  nearly  sonic 
slewing  velocities  do  not  occur  near  the  beam*  s point  of  origin. 

For  convenience,  Table  I lists  the  symbols  used  in  this  paper,  although 
the  variables  are  also  defined  when  they  are  first  used  in  the  text. 


- 6 - 


TABLE  I.  DEFINITION  OF  SYMBOI5 


Hydrcxlynamic  Variables 


p mass  density;  is  aunbient,  unperturbed  density 

r density  perturbation  divided  by 

R(k,t)  Fourier  transform  of  r(x,t) 

R(k,s)  Laplace  transform  of  R(k,t) 

p fluid  pressxire 

e si>ecific  internal  energy 

Y ratio  of  specific  heats 

c speed  of  sound 

V fluid  velocity 


Independent  Variables 


spatial  coordinate  measured  from  beam  center  at  t=0 

parameter  appearing  in  Laplace  transform 

spatial  wavenumber  appearing  in  Fourier  transform 


Properties  of  the  Beam 


X (t) 

o' 


position  of  beam  center 
beam  radius 


spatial  coordinate  measured  from  beam  center 
instantaneous  velocity  of  beam  with  respect  to  fluid 
Mach  number  of  beam,  V/c 


total  beam  power 


^ i.- 


TABLE  I.  (Continued) 


I(x,t)  beam  ir radiance 

0 normalized  irradiance 

E(k,t)  Fourier  1 ransi’orm  of  0 

^(kjr;)  I/iplace  trancform  of  F 


Heating  Rate 

a atmospheric  at.tenuation  coefficient 

Q rate  of  energy  deposition  per  unit  volume  of  fluid 

q normalized  (dimensionless)  rate  of  energy  deposition 


Time  Scales 


hydrodynamic  time,  a/c 

time  for  reaching  steady-state  perturbation 


T 


beam  acceleration  time-scale 


2.1  Basic  Equations 

The  fundamental  equations  describing  the  interaction  of  the  laser  beam 
with  the  atmosphere  are  (cf.  Ref.  7) 


V • ? .0 

(1) 

.0 

dt  p 

(2) 

p|^[|v^  + c]+V*  = Q, 

(3) 

vhich  represent  conseirvation  of  mass,  momentum  and  energy,  respectively.  Hie 
fluid  density,  velocity,  £ind  pressure  are  denoted  by  p,  i?,  and  p;  Q(i?^t)  is 
the  heating  rate,  and  is  given  by 

Q = al(x,t),  (4) 

idiere  a is  the  atmosi*ieric  absorption  coefficient  emd  I is  the  laser  irradiance. 
The  specific  Internal  energy  of  the  gas,  e>  is  given  by 


e 


£ 

P ’ 


(5) 


where  y is  the  ratio  of  specific  heats.  Using  the  above  equations,  equation  (3) 
is  easily  put  into  the  form 

f -vff  - (V-I)OI  (6) 
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LwearAaad  lifayaU<?ns 

The  time  derivatives  appearing  in  these  equations  are  Lagrangian 
derivatives,  i.e.,  in  operator  form 


dt 


y . 


(7) 


For  small -ajnplitude  perturbations  of  an  initially  uniform  fluid  at 
rest,  to  first  order  the  "v  • y terms  may  be  neglected  and  d/dt  =“  d/dt 
(Ref.  8).  Considering  now  only  fluid  flow  in  one  spatial  dimension 
denoted  by  the  coordinate  x,  the  linearized  forms  of  equations  (l),  (2)  and 
(6)  become 


dt 


dv 
*^o  dx 


(8) 


P 


^ ^ = 0 
o dt  dx 


(9) 


d£ 

9t 


7:1 


(10) 


where  second  and  higher  order  terms  have  been  omitted;  denoting  the  unperturbed 
pressiire  and  density  by  p^  and  respectively,  the  speed  of  sound  in  the  un- 
disturbed gas  is  given  by 

c = (\  • (11) 

In  equations  (8)  - (lo),  euad  throughout  the  following  development,  dependent 
variables  without  subscripts  denote  small  perturbations;  for  example,  the  total 
density  is  given  by  + p,  where  jp]  « p^. 
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Vftien  the  beam  moves  at  constant  velocity  with  respect  to  the  unperturbed 


air,  it  is  convenient  to  introduce  a new  spatial  variable, 

I = X - x^(t),  (12) 

where  x^(t)  is  the  position  of  bean  center  at  time  t.  For  definiteness,  we 
take  the  origin  of  spatial  coordinates  to  be  the  beam  center  at  t=0  (when  the 
beam  is  turned  on)  and  assume  the  beam  moves  toward  positive  x.  Dien,  with 
= Vt,  the  derivatives  with  respect  to  t in  equations  (8)  - (lO)  are  trans- 
formed as  follows*: 


St  St  " ^ 5f  • 


(13) 


The  (5,t)  coordinates  are  particxilaxly  convenient  for  describing  the  pertur- 
bation within  the  beam. 


If  V is  not  constsmt,  the  transformation  (12)  is  not  useful  for  two 
reasons.  First,  the  resulting  partial  differential  equations  contedn  terms 
in  which  derivatives  are  multiplied  by  the  non-constant  coefficient,  V(t). 
Second,  because  (?,t)  is  an  accelerated  reference  frame,  the  quantity  dV/dt 
would  appear  on  the  left-hand  side  of  the  momentum  equation;  since  it  does  not 
possess  a (spatial)  Fourier  treuisform,  it  would  be  necesseoy  to  take  the 
gradient  of  the  mcxnentum  equation  in  order  to  apply  this  transform.  Our 
procedure  is  to  solve  for  the  density  perturbation  in  (x,t)  coordinates,  then 
to  use  equation  (12)  to  evaluate  the  solution  in  the  region  of  interest  (l.e. , 
for  relatively  small  values  of  5). 

Note  that  the  wind  blows  with  velocity  -V  with  respect  to  the  bean. 
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2. Inhoitiogcneoua  Wave  Equation  for  the  Density 

From  equations  (fl)  - (lO),  it  is  easy  to  derive  a single  equation 
for  the  density: 


2 

c 


bi. 


(v-1) 


— 5 (O^)* 

ax 


Integrating  with  respect  to  time,  we  obteiln: 


(y-1) 


t 

2--  r al(x,f)df, 
dx  o 


(14) 


(15) 


using  the  conditions 

P = I*  = = 0 at  t = 0,  - ® < x < “ (l6) 

at^ 

(Ref.  lO).  Equation  (l5)  clearly  reveals  that  the  density  perturbation 

« 

is  a superposition  of  waves,  with  the  laser  heating  acting  as  driving  force. 
Physically,  energy  absorbed  from  the  beam  produces  local  expansion  which 
generates  waves  traveling  outward  at  the  speed  of  sound.  At  time  t^  after 
the  beam  is  turned  on,  p(x,t  ) is  non-zero  over  the  interval  -ct  s x s ct  . 
However,  the  solution  is  not  a superposition  of  sinple  (adiabatic)  sound 
waves  because  of  the  input  of  energy  to  the  gas  within  the  beam. 


* Our  solution  to  this  equation  presumes  that  l(x,t)  is  a known  function;  this 
is  a necessary  first  step  in  calculating  the  actual  intensity  of  a laser 
beam  propagating  throu^  perturbed  air. 


- 12 


Before  solving  equation  (15),  it  is  convenient  to  introduce  dlfflenslon- 
less  variables  — physical  quantities  are  measured  in  units  which  axe  paxti- 
cularly  relevant  to  the  problem,  and  continual  repetition  of  c , (y-l),  OL, 
etc.  is  avoided.  Introducing  the  units  defined  in  Table  II,  equation  (15) 
becomes 

q^  f 0(x,f)df,  (17) 

' dt  3x  / dx  o 

where,  for  example,  r(x,t)  denotes  the  actual  density  perturbation  divided 
by  the  ambient,  unperturbed  density,  p^. 
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TABLE  II.  DIMENSIOMLESS  UNITS 


QUANTITY 

UNIT  OF  MEASURE 

VALUE  AT  10  KM  ALTITUDE 

r,  Density  Perturbation 

p . Ambient 
o 

-k  -r 

~ 4 X 10  g cm 

M,  Velocity 

c,  Sound  Speed 

c ~ 300  m sec”^ 

X,  § Distance 

a.  Beam  Radius 

Assume  a = 0. ^ m 

t,  Time 

tR  = 

_3 

tjj  ~ 2 X 10  sec 

Beam  Irradlance 

CVJ 

II 

o 

M 

I^  =“  1.27  P 

q,  Energy  Deposition 

(Y-l)oI^tjj 

q _ g 

Po'' 

q ~ 2 X lO"^  aP 

= Energy  absorbed 
in  time  tjj 

^ Ambient  enthalpy 

for  a in  km“^ 
P in  kw 

r: . 4 Intcfjal  Tranarormations  of  the  Wave  Equation 


To  cari^  out  the  solution  of  equation  (l?),  we  apply  Fourier  and  Laplace 
transforms;  to  define  notation^  the  Fourier  transform  of  a function  g(x,t)  Is 
given  by 

C0 

G(k,t)  3 I*  g(x,t)dx,  (l8) 

•00 

and  the  Laplace  transform  of  G(k,t)  is  denoted  by 
00 

G(k,s)  3 f e"®^  G(k,t)dt,  (19) 

o 

Applying  these  transforms  to  equation  (I7)  results  in 

[s^  + k^j  R(k,s)  = -qk^  l'(k,s)/s.  (20) 

This  equation  is  valid  if  the  functions  r and  ^ (idsose  transforms  are  and 
respectively)  satisfy  certain  continuity  and  boundedness  conditions  (cf. 

Ref.  2).  To  obtain  a solution  we  assume  that  the  boTindary  conditions 
[eq.  (16)]  and  the  nature  of  the  physical  problem  prevent  the  appearance 

of  subtle  mathematical  pathology. 

Formally,  the  solution  for  p(x,t)  is  given  by  the  inverse  trsuisform  of 

R(k,,) . , (21) 

s(s  + k ) 

where  F is  presumed  known.  For  a beam  with  constant  intensity  profile  and 
>diicb  moves  such  that  its  center  at  time  t is  located  at  x^(t), 

fj(x,t)  = 0(x-x^). 
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:v,- 


(22) 


Using  the  translation  property  of  the  Fourier  transform  we  have 


F(k,t) 


-ikx  (t)  " 

o'  |« 


e 0(x)  dx. 


(23) 


where  ^(x)  gives  the  spatial  distribution  of  intensity  as  a function  of 
distance  from  beam  center.  A Gaussian  beam,  for  example,  is  described  by 

jZl(x)  = exp(-x  ) (24) 

(recall  that  distances  are  measiu-ed  in  units  of  the  beam  radius,  a).  To 
obtain  1'(k,s)  from  equation  (23),  it  is  necessary  to  take  the  Laplace  trans- 
form of  exp[-ikx^(t)] . 
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2 . 'j  Solution  for  tho  Penalty  Perturbation 

To  seek  a more  fruitful  approach,  we  go  back  to  equation  (2l),  expeind- 
ing  the  denominator  in  partial  fractions: 


K(k,s) 


2 


_1 

s-ik 


+ 


1 

s+ik 


(25) 


Hie  convolution  theorem  for  the  Laplace  transform  enables  us  to  write  the 
Fourier  transform  of  r(x,t)  as  a convolution: 
t 

R(k,t)  = q J F(k,f)  G(k,t-f)  dt',  (26) 

o 

vdiere  G(k,t)  is  the  inverse  Laplace  transform  of 


G(k,s)  = 5 j ’ s I ' 

explicitly, 

1 / ikt  . „-ikt\  T 

G(k,t)  = ^(e  +e  )-l 


(27) 


(28) 


and 


R(k,t)  = I q ,f  F(k,f)  dt* 


(29) 


Because  of  the  translation  property  of  Fourier  transforms,  for  anj-  given  values 
of  t and  t* , F(k,t’)  e*  is  the  Fourier  transform  of  ± (t-t*),  t*  ]. 

Therefore,  taking  the  inverse  Fourier  transform  of  equation  (29),  we  obtain 

r(x,t)  = i q j*  + (t-t*  ),  t*  ] + ^[x  - (t-t*  ),  t’  ]l  dt’  - 

o * ^ 

t 

q J 0(x,t*  ) dt*  , ( 50) 

o 

which  is  the  formal  solution  for  the  density  perturbation. 
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The  meaning  of  equation  (50)  becomes  cleeirer  if  we  consider  a beam 
described  by  equation  (22);  then 

0[x  ± (t-f  ),  t*  ] = 0[x  i (t-f  ) - xjf  )J.  (31) 

The  expression  x ± (t-f)  represents  the  position  at  time  f < t of  a sonic 
impulse  which  will  pass  the  point  x at  time  t;  the  minus  sign  corresponds 
to  an  impulse  moving  in  the  positive  x direction  (a  "forward-facing  charac- 
teristic") and  the  plus  sign,  to  a wave  moving  toward  negative  x ("backward 
facing  characteristic").  The  argument  of  0 in  equation  (31)  is  therefore  the 
distance  (at  time  f ) between  the  beam  center  and  the  cheu'acteri sties  passing 
through  the  point  (x,t),  and  the  first  two  terms  on  the  right-hand  side  of 
equation  (30)  represent  the  effect  at  (x,t)  due  to  energy  which  was  absorbed 
at  points  on  the  characteristics  during  the  Interval  0 to  t.  Figure  1 is 
a space-time  diagram  which  Illustrates  the  geometry  involved.  Ihe  dashed 
curves  on  either  side  of  x^(t* ) represent  the  spatial  extent  jf  the  beam, 
within  which  power  absorption  is  greatest.  As  shown,  the  point  (x,t)  has 
been  within  the  beam  for  several  hydrodynamic  times  (since  f =“  6.5)  and 
the  resultant  absorbed  energy  contributes  an  expansion  of  the  air,  given  by 
the  last  term  of  equation  (50).  In  addition,  air  with  4 ^ x s.  6 has  been 
within  the  beam  at  times  such  that  effects  of  the  perturbation  can  propagate 
to  (x,t).  Adding  energy  to  a small  volume  of  air  produces  local  expeuislon; 
the  "snowplow"  effect  on  neighboring  volume  elements  produces  a positive  denr.iiy 
perturbation  outside  the  originally  heated  volume.  3he  net  perturbation , as 
given  by  equation  (30)  is  therefore  a superposition  of  these  two  effects. 
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2.6  The  Constant  Velocity  Bead 


As  a special  case,  we  can  evaluate  equation  (JO)  for  a beam  moving  at 
constant  velocity  (or,  equivaloitly , a stationary  beam  and  a constant  trans- 
verse wind).  Put 

x^(t)  = Mt,  (52) 

where  M Is  the  mach  number  for  the  tmperturbed  flow.  Using  equation  (Jl), 
i.e. , assuming  the  intensity  of  the  beam  is  constant,  the  integrals  for 
r(x,t)  are  of  the  form 

t 

J 0 [x  ± t - (M±  l)t'3dt«.  (33) 

o 

Defining  a variable  y=xit-(M±  l)t’  , (33)  can  be  written  as 
x-Mt 

-whi  f 

x±t 

For  a Gaussian  beam  (cf,  eq.  [24]),  the  integral  is 


/tt  erf(x±t)  - erf  (x-Mt) 
T"  Mtl 


(35) 


(for  properties  of  the  error  function,  erf(z),  see,  e.g..  Ref.  l).  We 
are  primarily  interested  in  the  density  profile  within  the  beam  itself, 
so  introducing  the  displacement  from  beam  center  (§  is  defined  in  Eq.  (12)), 
the  solution  becomes 

t)  - ±(M-H)tI  1 erf[§  ±(M-l)tJ 

ns,t;  - 2 |2  M + 1 2 M-1 


erfCS  ^ 
M 


er 


M(  Mtl) (M-1) 


(56) 
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ITiis  fuiu'lion  i;:.  :'.ketcht,’<J  in  Figure  2 -i'.  a time  ita-  wliich  i.iic  ()roJ  Lii-  wi ; iii  n 


t.he  betun  ia  iieaJ'  ita  s1.ead,y -rotate  shape.  iiie  poird  •- > (ii  ijirun  <■  nii.-f)  mov< 

to  t he  rigtit  at  speed  M,  ieavinf'  behind  a low-d(  ri:;l  ty  rer;i,,ri  „/bi'-h 


ime  int  erval,  ttie  true  ste-ad;/-sl.at«  den:u  !y  is  reach«  d as 


lim 


M(i-M)(l+M) 


a.a-l  u riraj 

i;: 

1 ( i • , 

’ ti'.'  besun 

a',-  a;'t.'-r 

a I'ini'  e 

as  ■ 

for  ■ 

2 _ 

(f7) 

S'l  Sisaoy 

i nf  ini  i.y 

'•annf,  se 

assumed  to  be  unperturbed.  Fr^m  •.•puaviou  ( ly'f)  the  stead,v-s' a'  * aonuitir.ns  at 
si)atiaL  infinity  are 


/ ^ /tt  q 

= 2(W  • 


The:,!  non-zero  values  are  due  to  the  fact  ttiat  in  i,ne-dimensii;nai  flow  without 
disi'.iputioii,  the  amplitude  of  a wave  docs  not.  dec-rease-  with  distance  from  tlie 
source;  for  flow  in  two  <tr  t.liree  dimensions  we  would  expi'ct  lo  have  o-ro  per- 


turbation at  infinit.y. 


X 


J'or  ;5upc:rsonic  i'low  Nf'  1,  and  ojuation  (jXO  itoad^y-i^tal' 


limit. 


r(?, 


Ar  q 1 - 

m(m+i)(m-i) 


(A) 


Ttie  air  :iri  front  of  Uic  boaiti  {.nxiy  Z)  unporturbod  L';uauG('  jiroGGuro  waves 

created  liy  ttie  be.ai.'i  propa<;at<;  more  slowly  ttian  the  Learr,  itself. 


Equation  (yo)  also  exprt.'sses  the  solution  for  the  sp'.el.al  case  M = 1; 
proceccllnf;  as  bi,!forc,  one  easily  obtains  tlie  time-dfqjondenl.  .jc.lution 


(t  e^(-  f ) + I 
( /tj  '2 


rf(§-K_'t) 


frf(§) 


I 

) • 


In  tiiis  case  there  is  no  steady-state.  For  t > 1 t)ie  first  term  of  equa- 

tion (^J)  dominates  near  beam  center;  it  describes  a Gaussian  density  profile 
with  maximum  at  beam  c;c'nter.  The  rrafyiitude  of  tlie  perturbation  frows  linearly 
in  time.  Hie  present  solution  remains  valid  as  lonf  as  |r(P,t)|  « 1,  wtiich 
in  turn  itr^jlies  t l/lOq.  Gince  expected  values  of  q :.a'o  muc)i  smaller  than 
unity  (see  'Table  II),  (,'quation  (*40)  describes  the  sonic  case  for  significant 
period.s  of  time. 


At  low  wind  speeds,  the  approach  to  a steady-state  witliin  the  beajn  is 
propxjrtionai.  to  the  wind  transit  time,  wliich  varies  as  M ^ (cf.  Ref.  11). 

At  )iigher  speeds,  for  M ~ 1,  sufficient  time  is  required  for  the  forward- 
moving  wave  to  proceed  out  of  the  beam.  Tiiis  wave  moves  realtive  to  the 


beam  at  mach  nvUDbcr  1-M,  so  the  approach  to  steady "Stati,;  is  longer  as  M 
increases.  Quantitatively,  the  profile  given  by  equation  is  within 

1 percent  of  the  steady-state  for  -2  ^ | s 2 when  t s with 

At  max  {4  M 4(1-M) 
eq 

As  can  be  seen  from  Figure  3>  the  steady-state  density  change  across 
the  beam  (from  ? * -2  to  ? = +2)  becomes  large  as  M approaches  0 or  1, 
although  the  time  required  to  reach  the  steady-state  (At^^)  also  increase^. 
Tiirtc -dependent  density  profiles  given  by  equation  (3^)  arc  in  complete  agree- 
ment vrith  Munn's^'^^  results  for  the  constant  velocity  beam. 
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P. . 7 Tlie  laobarle  Approximtion 

l''ur  F.uf ficit  lit  iy  small  M,  1,h<  density  perturbat  ioii  is  produc'  d primarily 
by  expaiisiun  at  nearly  ' ons i.arit.  pressure  and  the  t-ulf  nV  ijutraiiij'  waves  is  small. 
The  approxinval  e r.olui  ion  is.  obi.ained  from  equat  ion  ( lo)  wit  ti  rb’  [iresstire 
perturbation  set  identically  t.o  zero,  and  is  dveri  by  itie  third  term  on  the 
right-hand  ride  of  eqiiat  ion  (iXl).  Hie  isobaric  rolut  ien  ' or  (iaussi'iri  beam 
is 

r - yjUi  erf(F)  - erf(g  + Mt)  , ^ > 

The  corresponding  steady-state  density  change  across  the  beam  is  shown  in  Figure 
3 by  the  dashed  curve.  The  exact  solution,  equation  (>6),  contains  terms  whicn 
are  smaller  than  those  in  equation  (42)  by  a factoi-  of  approx  i irately  M, 
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2.8 


C'ompariGon  of  Constant  Velocity  axid  Acceloratcd  Beams 


In  the  general  case,  vdien  the  beam  moves  with  arbitrary  velocity,  the 
Integrals  appeai’ing  in  equation  ( JK))  must  be  numerically  evaluated.  Note  that 
the  beam  intensity  may  also  vary  arbitrarily,  so  pulsed  beams  can  be  treated 
as  well  as  CW. 

Figure  4 shows  the  density  perturbation  for  a Gaussian  beam  moving 
at  K - 0.1;  when  t = 40,  the  profile  has  essentially  achieved  a steady-state 
for  -2  s F s 2,  vAiere  the  beam  Is  ra<jst  intense.  Figure  5 shows  the  profile 
for  a beam  vrtiich  initially  moves  at  M = 0.1  and  is  accelerated  at  the  rate 
dM/dt  = lu  ^ (in  dimensionless  units);  for  conditions  at  10  km  altitude 
(cf.  Table  II)  this  corresponds  to  an  acceleration  of  I50  m/sec  . At  early 
times  the  perturbation  does  not  greatly  differ  from  the  constant-velocity 
case,  but  by  t = 100  the  departure  is  obvious.  Between  t 100  and  t =“  "JOG 
the  perturbation  closely  approaches  the  steady-state  profile  corresponding 
to  the  instantaneous  velocity.  The  later  development  is  shown  in  Figure  6. 

As  M exceeds  0.8,  there  develops  a density  maximum  near  § = 1 and,  especially 
for  § > 0,  the  profile  departs  more  and  more  from  the  corresponding  steady-state 
in  the  sense  of  having  smaller  perturbation.  As  a measure  of  the  perturbation 
at  eeich  instant,  we  use  the  dlfferenceiln  density  frean  maximum  to  minimum  with- 
In  -2  5 '•  2.  The  solid  curve  in  Figure  7 represents  the  magnitude  of  the 

steady-state  perturbation  for  constant -velocity  beams  as  a function  of  mach 
number  (from  Figure  3).  Uie  points  connected  by  a dashed  curve  represent  the 
profiles  shown  in  Figures  5 “id  across  the  top  of  Figure  7 is  the  time  scale 
for  the  accelerated  beam.  When  the  beam  is  turned  on  at  t = 0,  M = 0.1  and 
the  density  perturbation  Is  zero.  At  this  speed,  a steady-state  would  be 


- ^ - 


t^iguro  b.  Time -dependent  density  profiler  for  an  accelerated  beam; 

initially  M = 0.1,  dM/dt  = 0.001.  Stor  t a ^300,  nee  Fig.  6 


1 


TIME 


l''i(5ur«'  Y.  l\'n;iL(.,y  chunj'/-  across  ai'<  > lATatid  l)oam  (point/3  and  da.'.ti'  d curv"; 
ijomparod  with  .'it«.'ady-3tal'  condition.!  for  nonjiant  vtlf'clty 
bi-am.)  (.Jolid  curve) . 
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1 


fcacht'd  when  1 but  by  i ha*  time  ’he  beam's  mach  number  has  increased  to 

O.lti,  Tile  new  ct'eed  Irij  lie:  i :'Ji  irt  ei-  •■■jui  libmt  ion  time  and  a smaller  steady-state 

perturbation;  ’o^’ethei"  wl’h  the  i'Jt.i'e  e.jui  1 ibrat  i' )H  time,  this  produces  a 
slightly  l.atve'-  ner’aflia’i  't  'i.ai.  w.  i.:,i  le  • tn  -atte  f- f 'he  instantaneous 
steady-state.  As  M.e  i eai:  Vf-s  : a:  ’•  i , ■ :.<  ,;ui  J it  :‘<i' lu.  lime  decreases 

and  the  ner*  irt,.,- j , p.  : ea-  r.  <» , r;t  efidy -s  i.at.e.  Later, 

when  M ► e.  , '.he  e,-ui  j 1 1 .-a  • : : . . ■ ,.f  magnitude  of  the  pertur- 

bation lags  betilnd  'he  ■ • e.au,. - 

Tile  pt'S’ urhativ.n  l ei  i >i.- -e i .t  i' • ,-;ur.  is  s.tiown  in  Figure  B.  Initially, 

M = 0.1;  with  dM/dt  -Lo  , 'tie  t'etur.  is.  fa'i'iiar;/  at  t = 100  and  moves  back 
toward  its  initial  position  at  later  'imes.  Wiien  t ^ 100,  the  amplitude  of 
the  perturbation  is  less  than  Uie  steady-state  value,  but  after  t = 100,  the 
perturbation  considerably  exceeds  the  steady-state  value.  At  t = I50,  with 
M = -O.O5,  the  steady-state  density  change  across  the  beam  would  be  approxi- 
mately Ar  3t<q,  while  the  actual  Ar  is  nearly  three  times  as  large. 

Defining  a time-scale  for  acceleration  by 

t=1m/a1,  (4>) 

vdiere  A = dM/dt  (in  dimensionless  units),  it  is  possible  to  estimate  those 
velocities  and  accelerations  for  which  the  density  profile  will  differ  signi- 
ficantly from  the  constant-velocity  solution.  When  t is  smaller  thsui 
the  beam  velocity  changes  more  quickly  than  the  density  approaches  steady-state. 
Using  equation  (4l),  this  condition  becomes 

1a|  > min  {M^/4,  (l-M)^/4},  (44) 

from  which  we  see  that  small  accelerations  have  the  greatest  effect  for  M “ 0 
or  M 1. 
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5.0 


TWO-DIMENSIONAL  THEHMAL  BLOOMING  FOR  ARinTRARY  TRANSVERSE  WIND 


It  has  proven  convenient  to  part  from  traditional  formulations  of 

the  laser  heating  hydrodynamics  which  assame  a simple  relationship  between 

Eulerian  (beam  centered)  and  Lagrangian  (parcel  of  air  centcrc^d)  coordinate 

systems.  The  treatment  given  here  does  not  assume  t)i(,'  transverse  I'iow  to  be 

(l  2) 

constant  in  magnitude  or  direction.  ' ' 


In  Eulerian  coordinates,  at  a fixed  distance  7.  along  the  beam  the  wind 
flowing  in  the  x,y  plane  (perpendicular  to  the  ::-axis)  is  assumed  to  be  a 
known  function  of  time: 

w(z,t)  = X (z,t)  + y Wy  (z,t)  (1) 

and  in  general  w is  composed  of  coD5>onents  due  to  natural  wind,  motion  of 
the  laser  platform,  and  non-zero  slowing  of  the  beam.  The  La^'rangian  formu- 
lation consists  of  iixprossing  the  e^^uations  of  hydrodynamics  in  the  coordinate 
frame  which  is  at  rest  in  the  unperturbed  air.  In  this  frame  the  beam  moves 
with  average  velocity 

v^(z,t)  ■ - w (/..L) 

In  the  discussion  below,  note  that  all  equations  are  written  in 
physical  variables;  dimensionless  variables  are  not  introduced. 


- j> 


MatlicrriatJ  cal  JoJutLon 


In  Lo/'iranp.ian  <-oordirualc3 , Jn  wti.irh  Ui'-  uri]ji  J't.uxlicd  !<:•  ii\w.  i,;  il 
nc'.t,  ciaaJl  amjilitudn  varjations  in  ]'reJr,ur' ■,  di  n.'.iiy  and  vi  ii  ^iLy  (d'url'.'d 
Ly  j),  p and  v,  reij>ectiveiy)  natinfy  the  iineax’i /.ed  equation.;; 


^ +0  V • V = 0 

■ )t 


t = (y  - 1)  a I (7,  t) 

ot  p ot 


(i) 

(4) 


(!^) 


Tilo  ambient  })ros3\xre  and  density  are  and  p^,  Y in  the  ratio  of  sjjccific 
heats  and  a is  the  atmospheric  absorption  coefficient.  Ihe  beam  irradiance 
I includes  time  dependence  vriiich  is  due  to  motion  of  the  beam  relative  to 
the  unperturbed  medium. 


Intetyrating  equation  (‘y)  with  respect  to  time  gives 


pQ'^  ^ t)  _ (v  -•  j) 


a 


YP 


YP. 


J I (7,  t ) (it 

^ f > / • 


(-:■) 


On  the  ri^ht-hand  side,  the  first  term  represents  an  acoustic  (adiabatic) 
disturbance,  wtiile  tin;  second  term  is  an  isoljaric  expansion  whicii  dominates 
when  the  boEim  moves  slowly.  Uie  pressure  variation  p can  itself  be  expressed 
as  an  integral  of  the  irradiance  multijilied  by  an  appro])riate  Greenh  Itincticri 
wh(;reupoti  substitution  Into  equation  (<■)  gives  tin'  solution  for  p In  term.; 
of  known  lYinctions. 
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Into  equation  (4)  we  introduce  the  vector  potential  i)f(r,  t),  which 
is  related  to  the  velocity  by 

V = - (l 

It  is  then  easily  shown  that 
-E_  _ d!L  , 


where  c = (yp^/Pq)  ^ is  the  sovind  ;^eed.  As  follows  from  equations  (3), 

(4)  and  (?),  ijr  satisfies  an  inhomogeneous  wave  equation  with  source  term 

(y  - l)  Ql/p^c  , and  it  can  be  written  in  terms  of  a Green's  function 

G(r  ,t  ; r,t)  which  gives  the  contribution  at  (r,t)  due  to  energy  deposited 
o o 

at  (r  ,t  ).  Bie  time  derivative  of  ♦ is  then 
o o' 


J I (r  » t ) ^ 

< " oj  O 'ci'r.'rt+. 


o’  o’  dt 


Morse  and  Feshbach^^^  give  the  appropriate  Green's  function  for  the  two- 
dimensional  problem;  i.e.,  for  ^atial  variations  in  the  (x,  y)  plane  normal 
to  the  beeun  and  neglecting  slow  variations  along  the  beam  (parallel  to  the 
z axis).  The  spatial  integral  in  equation  (9)  then  reduces  to  a .surface 
integral  which  can  be  expressed  in  polar  coordinates  (R,  6)  centered  at  the 
point  (x,  y) : 

R H (x  - x^)^  + (y  - y^)^ 

0 = arc  tan  l(y  - y^)/(x  - x^)]  • ^ 


The  result  can  be  put  into  the  form 


p c 

o 
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1 


(11) 


where 


F(e,  ; X,  y,  t)  - i(R^,  0.  t^) 


R 


- I 


max 


R dR  

° 1/ 


ifH,  e,  i ) - i{R  , e,  t ) 

__  O m/iX  O 

57^ 


\{t  - i r - R'  /r 
I ()  ' 


.r 


(12) 

and  “ c(t  ~t^)*  By  choosing  the  coordinates  (x,  y)  to  lie  within  the 

beam  at  time  t,  the  solution  is  evaluated  only  at  points  of  interest  for 
calculating  thermal  blooming. 


1 


J.2  Compajiaon  with  the  ftie -Dimensional  Solution 

The  Green's  function  formalism  for  expressing  laser  induced  density 
changes  makes  it  somewhat  difficult  to  develop  an  intuitive  understanding  of 
the  underlying  physics.  In  order  to  deomonstrate  use  of  the  formalism,  and  to 
establish  contact  with  a known  and  relatively  easily  solved  case,  we  consider 
the  one -dimensional,  constant  velocity  case  discussed  in  Section  2.  For  this 
case  equations  (6)  and  (9)  become,  with  Q *=  cOi, 

. t . oG  I 

p(x,  t)  = - J dt'  q(x,.  t')  - J dx'  q(x'.  L)  ^[(13) 

c*'  o ' c 

and  the  Green's  function,  given  by  Morse  and  Feshbach,  is 


wIk  T'-  >)l‘  lii-.'uii  (•<  tiLcr  at  f.irtic  t..  l''or  th'  <jajc 

or  (■( l.r.-ui.:vi 'Imo  v/.IikI  wc  ni.a,y  put 

= vt  • (.1  ;) 

Note  that  tii('  wind  riK'Vor.  at  .spei’d  w = - v with  I to  the  t'rain  and  th  ■ 

intcfirai.'.  a]ip<‘ariar.  In  ouiuation  (ih)  arc  of  the  I’ortii 

t 

J dt'  Q,(x  i c[t  - t']  - vt')  (id) 

o 

which  can  be  ninplified  by  the  substitution 

y = X ± c[t  - t']  = X t ct  - (v  ± c)t'  ' (20) 

Tlicn  (19)  oncotnes 

• 7TT  i*  , , ^ (21) 

xtet 

To  express  this  result  in  coordinates  which  move  with  the  beam,  mai'.o  the 
transformation 


5 = X - vt 


(22) 


F is  the  distance  from  beam  center.  Equation  (2l)  ihen  become: 


/±  cl.  , ^(y) 

?+(v»c)t 


(25) 


Written  out  exjjlicitiy  in  terms  of  the  wind  speed  w,  equation  (16)  now 
becomes 


Q(y)dy  I 

(c+w)t  '^-'"^^.(c-w)t 


Q(y)dy 


(24) 
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when  |w|  «<  c the  2nd  and  3rd  terms  are  negligible  with  respect  to  the  1st 


term  stnd  we  have 


p(5.t) 


(Y.  -..1)  g.  r i;(y)  ^ 

c w c-wt 


(25) 


which  is  the  familiar  form  of  the  solution  for  the  isobaric  approximation. 


- 41  - 


50 


Density  Variations  within  an  Accelerated  Beam 


Using  the  equations  developed  in  Section  5.1  we  have  calculated  the 
thermal  blooming  density  perturbation  for  an  infinite  Gaussian  irradiance 
profile.  Tb  illustrate  the  effects  of  propagation  in  the  presence  of  stagna- 
tion zones,  as  in  Section  2.8,  we  chose  the  time -dependence  of  the  beam 
motion  to  be 

V - - at  (26) 

with  initial  velocity  v^  = 0.1  c and  constant  acceleration 
“3  2 2 *1  2 

a = 10  c /o>  ~ 10  uj  m/sec  , where  cn  is  the  beam  radius  (in  meters).  The 
flow  velocity  relative  to  the  beam  is  w = - v;  at  a time  At  after  beam  turn- 
on, the  flow  velocity  is  momentarily  zero,  vrtiere 

At  a 0.5  a)  sec  . (27) 

Figure  1 shows  density  contours  for  t s 1.5  At.  During  this  interval 
the  density  never  reaches  a steady  state,  nor  can  the  profiles  be  approximated 
by  suitably  chosen  steady  state  curves. 


Figure  1. 


Density  Contours  In  an  Accelerated  Bmud.  At  <,ach 
time  t are  shown  curves  of  constant  |dp/Pol  meacurca 
in  units  of  (Y-l)0flElD/poC^.  The  spatial  coordinate:-. 
(5/Tl)  are  co-moving  with  the  beam;  dlstences  are 
measxired  in  beam  radii,  cd.  The  flow  velocity  has  the 
same  magnitude  at  t « 0.5  dt  and  t » 1.5  dt,  but  the 
density  profiles  are  sig^ficantly  different  both  in 
shape  and  magnitude. 
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IRRADIANCE  PROFILEi;  FOR  NON-COPLANAR  CLEW  AND  WIND 


Ij.O 

To  protLiot  hho  irradiance  within  a la;;or  beam  propa(^tinp  through  tno 
atmocriilii'ri:  ri’ciuirca  a wav<  o]>tjnn  comi)ut(.-r  code  which  includes  an  evaluation  of 
the  laocr-Liiduced  )iydro dynamic;;  which  war.  dincucacd  in  earlier  cectionc. 

Dr.  C.  h.  Ho,;":'-  of  the  Air  Force  Weapon;;  lyiboratory  in  Albupunrciue  kindly  tnade 
.avai.l.ahJ. ■ to  u.;  hiu  wave  optic;;  code,  whiich  wc  modified  to  carry  out  the  com- 
putation;; diccuaned  below. 

Aa  a first  step  toward  in^lomontinc  the  mathematical  solution  derived 
in  Section  we  replaced  the  hydrodynamics  section  of  Dr.  Hogge's  code  with 
coding  wJiich  calcxolates  two-dimensional  density  perturbation  for  arbitrary 
magnitude  and  direction  of  the  transverse  flow  velocity,  asstiming  steady-state 
conditions  obtain.  Further  generalization  of  the  coding  to  deal  with  the  fully 
time -dependent  problem  is  now  relatively  straightforward.  But  before  making 
the  second  step,  we  made  some  tost  calculations  which  in  themselves  yielded 
useful  results. 

A number  of  runs  were  made  with  the  modified  propagation  code  for  cases 
in  which  the  natural  wind  (perhaps  duo  to  laser  platform  motion)  is  not  coplanar 
with  the  non-zero  slew  velocity.  From  the  results  we  have  found  that  a single 
scaling  is  possible  to  predict  the  peak  intensity  in  the  focal  plane  (under 
certain  conditions)  from  coplanar  data. 


- k-f  - 


4.1 


Gcenario  Definition 


Propagation  runs  were  performed  for  the  following  set  of  conditionr:: 

• Co-altitude  device  and  targf^t 

• Platform  induced  flow  in  the  lioriiiontal  plane 

• Slew  plane  oriented  at  arbitrary  angles  (0°  to  lTo°  ) witli 
reject  to  the  horizontal. 

A number  of  different 

• Power  levels 

• Beam  sizes 

• Kinetic  cooling  conditions 

• Slew  rates 

wore  considered.  In  particular  slew  rates  were  adjusted  so  that  the  l80° 
orientation  included  a case  which  would  have  had  a stagnation  zone,  although 
no  calculation  was  performed  for  this  case. 


At  0°  slew  and  wind  velocities  are  additive  and  at  l80°  they  subtract;  at 
intermediate  angles  the  vector  sum  has  a non-zero  conponent  perpendicular 
to  the  horizontal.  (This  conponent  is  a linear  function  of  distance  along 
the  beam. ) 
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4.2 


Scaling  Considerations 


A particular  subset  of  tlio  code  runs  were  found  amenable.'  to  scalirii^ 
for  determining  peak  intensity  levels.  The  conditions  apjiropriato  for  this 
are 

• Moderate  to  little  kinetic  cooling 

• Steady  state 

• Slew  angle  less  than  90° . 

Define  in  the  ustial  way 

= l(PEAK,  BLOOMED)/! (PEAK,  VACUUM,  ATTENUATED)  (l) 

KKlt 


where  I is  the  intensity.  Also  define  an  effective  transverse  flow  velocity 
as 

^EFF  “ I (2) 

o 


where 

v(z)  = V + nz 

' ' CD 


(5) 


and  R is  range  to  target. 


Here  is  natural  wind  and/or  platform  motion  induced  flow  and 
Q is  the  slew  rate.  The  vector  notation  formally  signifies  that  v^^  and  Q 
are  not  necessarily  coplanar. 


All  calculations  shown  sire  for  a beam  focused  on  a target  5 km  from 
the  transmitter  and  an  absorption  coefficient  of  0.1  km  the  platform 
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V''ir.ci1y  ij  V = ('7.0';  m/si.'c  (I'jO  nijih).  A.i  exainpie.;  cil'  t.hc  comjiu t-'iLj onal 
roault.s,  Figurf's  1-;’  ,:how  tlie  local -j)irun ■ irmtiLar'-e  criuLr.ar.;  Tor  G -■  0,  'i(f 
and  K'O.  ro.Tpectivt'ly . The  iKucimum  irradi.'uu'c  occur."  at  tl:o  pi'isit  ni.arkcd  hy 
tiie  CTO.:.;;  Die  contour.,  arc  for  the  'jOio  and  lO^  irracLLanco,  normlizcd 

to  t)ie  maximurii.  Bic  .;low  rate  is  0.0268  rad/sec,  and  the  cLLrection  of  t.he 
resultant  slc;w  velocity  cortponent  is  shown  in  each  figure.  This  anfptlar 
velocity  was  chosen  to  give  a slow  velocity  equal  in  magnitude  to  v at  the 
iiolf-way  point  along  the  beam.  Figures  4 and  5 show  the  mai_|iiitude  ol'  the  flow 
velocity  and  its  orientation  with  respect  to  v^  as  functions  of  distance  along 

the  beam.  Figure  4,  for  0 = 90°,  is  qualitatively  representative  of  B t,  j0° 

in  the  sense  that  the  transverse  flow  velocity  is  a monntonically  increasing 
function  of  distance  flxjm  the  transmitter.  In  contrast,  Figure  ‘j  shows  how 
the  flow  velocity  for  6 > 90°  (120°  in  this  case)  has  a minirman  between  the 
transmitter  and  target.  Bie  effective  v<;locity  (V  „,.)  is  d'jfiried  by  "Tuatiori 
(2),  and  is  the  average  along  the  beam  of  the  magnitude  of  the  flew  velocity. 

T>ie  centroid  of  the  focal  spot  is  displaced  from  the  optic  .axis  of 

the  transmitter  (the  aim  point  in  the  focal  plane).  Tliis  is  shown  in  Figure  6. 

Increasing  0 moves  the  beam  farther  from  the  aim  point  and  away  from  the 
horizontal  fax!;;. 

Figure  7 shows  as  a function  of  each  point  represents  a 

calculation  using  a modification  of  Dr.  Hogge's  wave  optics  code.  The  infinite 

Gaussian  has  a spot  size  of  0.25  m;  the  slew  angular  velocity  is 
■“2 

0^  = 2.682  x 10  rad/sec.  From  left  to  right,  the  points  represent 
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'iT  wi  -.d 


■^ure  plant  i rranci  i;inr'e  contours.  Wi-id  T.d  sltw  ('av<-  mcr.itcifs 


FLOW  DIRECTION  (degrees)  TRANSVERSE  FLOW  VELO^TTr' 


Figure  U,  Magnitude  and  direction  (measured  from  horizontal)  of  the  flow 
velocity  as  iXuictions  of  distance  along  the  b(;am,  for  I • VICP  . 
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r ill^d  r'^nr^ 


0 “ 15°,  50°  , *^5°  , 6cf  and  9<^  . !nie  truncated  Gaussian  has  a spot  size  of 
0.1768  m;  calculations  were  done  for  eund  for  = 0.5  The  coplanar 

cases  are  indicated  in  the  figure.  The  triangular  symbols  represent  non 
coplanar  cases  with  fi  ■ while  those  for  Ci  *=  Ci^  are  shown  as  square  symbols. 
Die  open  symbols  denote  cases  with  0 > 90° , while  filled  symbols  are  for 
0 ^ 90° . As  long  as  0 is  not  greater  than  90'' , there  is  a striking  lin ,‘ar 
relation  between  ^eff'  linear  dependence  resembles  the  usual 

scaling  laws  for  coplanar  wind  and  slew;  what  is  demonstrated  here  is  a linear 
scaling  of  vrtien  the  wind  and  slew  are  not  in  the  same  pleme,  and  the 

validity  of  the  scaling  to  6 > 90° . 

With  the  quantities  defined  at  the  beginning  of  the  section,  non 
coplanar  peak  intensities  may  be  obtained  as  follows; 

a)  For  a given  v , calculate  I , in  the  standard  way  for  two 

ui  rel 

coplanar  cases,  (v^,  with  ^ ^ 0. 

b)  Using  these  two  values  of  the  two  corresponding  values 

of  derive  the  straight  line  = A + B. 

c)  For  a non  coplansu:  case,  calculate  V „„  for  v(z)  = v + 5z 

ell  0) 

using  Equation  (2),  then  find  from  the  equation  derived 
in  Part  (b)  above. 

Diere  is  a caveat  attached  to  this  ^rtiich  must  be  strictly  observed,  namely 

a)  Meignitude  of  v(z)  must  be  monotone  increasing  (i.e.,  0 must  be 
no  greater  than  9CP  ) . 

More  con^jlicated  scalings  exist  for  other  cases  in  the  presence  of  strong 
kinetic  cooling  or  non  monotone  increasing  velocity  magnitude.;.  Scaling  for 
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these  cases  requires  more  parameters.  We  also  note  that  no  conclusions  should 
be  dravm  regarding  beam  shape  ajid  size  as  conpared  to  the  coplanar  case  since 
non  symmetrical  distortions  do  occur  and  the  relationship  between  peak  and 
average  intensity  changes. 


